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We provide a quantitative determination of the effective partonic kinematics for Higgs production 
in gluon fusion in terms of the collider energy at the LHC. We use the result to assess, as a function 
of the Higgs mass, whether the large mt approximation is adequate and whether Sudakov resum- 
mation is advantageous. We argue that our results hold to all perturbative orders. Based on it, we 
conclude that the full inclusion of finite top mass corrections is likely to be important for accurate 
phenomenology for a light Higgs with rnu ~ 125 GeV at the LHC with *J~s = 14 TeV. 



An accurate computation of the Higgs boson produc- 
tion cross-section [T] is essential in the search for this 
particle, which might be on the verge of being discovered 
at the LHC 0. If the Higgs is light, and in particu- 
lar in the region m# ~ 125 GeV, perhaps favored by 
LHC data, the dominant Higgs production mechanism is 
gluon fusion, which starts at leading order 0(a 2 ) through 
a (predominantly top) quark loop. Higher order correc- 
tions, which turn out to be quite large, are accordingly 
difficult to compute, and the full next-to-next-to-leading 
order (NNLO) result is known either in the large m t limit, 
at the fully differential level [3] , or as an expansion in in- 
verse powers of mt for the fully inclusive cross-section [3] . 
As mt — > oo the quark loop shrinks to a point (pointlike 
approximation) and the LO process becomes a tree-level 
process of an effective theory. 

At NLO, where the exact result is known, the large mt 
approximation turns out to work surprisingly well, even 
up to values of the Higgs mass at and above the top pair 
production threshold. This result can at least in part be 
understood based on the observation that the partonic 
cross-section is dominated by logarithmically enhanced 
contributions related to soft gluon radiation which are in- 
dependent of m t , so the pointlikc approximation becomes 
exact, up to an overall factor which starts at NNLO [5]. 
This soft dominance should take place when m# is raised 
at fixed s so the energy s of the partonic production sub- 
process approaches threshold s ~ m 2 H . 

Close enough to threshold it is advantageous to resum 
these logarithmically enhanced terms (threshold resum- 
mation), and it has indeed been observed that this re- 
summation significantly corrects and stabilizes the per- 
turbative result in regions in which the pointlike approxi- 
mation holds to satisfactory accuracy [5] . It is important 
to understand that this may happen even if the expan- 
sion in powers of the strong coupling a s behaves in a 
perturbative way, i.e. if the size of higher-order logarith- 



mically enhanced contributions decreases with the per- 
turbative order, so an actual all-order resummation is 
not really necessary. For this, it is sufficient that these 
enhanced contributions approximate the missing higher 
orders well enough that their inclusion actually improves 
the accuracy of the computation, and the desirability of 
resummation should thus be judged by the accuracy of 
the logarithmic approximation. 

Even so, this can only be possible if the center of 
mass energy of the partonic collision s is significantly 
lower than the hadronic one s, which at the LHC is very 
far from threshold. Because the gluon distribution is 
strongly peaked at small values of the momentum frac- 
tion of each hadron carried by the gluon itself, this is 
especially likely in a gluon fusion channel. A quantita- 
tive assessment of this effect is thus important in order 
to determine the accuracy of the pointlike limit, and also 
whether threshold resummation is advantageous. 

The necessary formalism was presented in Ref. [7] , and 
applied to Drell-Yan production. The cross-section for 
Higgs production is a function of the scale m 2 H and a 
scaling variable r = m 2 H /s, given by the convolution 

ajr, m 2 H ) = f 1 dz ^ /r ^ 2 \ ajz, a s (m 2 H )) 

T J T Z \z' H J Z 

of a partonic cross-section a , a s (m 2 H )j and a par- 
ton luminosity Jz? (| , mfj); a sum over relevant partonic 
subprocesses is understood. In our case at LO the only 
subprocess is gg — > H, so 

j^(z, M 2 )= [ ld ^9(y,n 2 )g(-,iA (2) 
Jz v \v J 

where g(z,^i 2 ) is the gluon parton distribution (PDF) in 
the proton. It is convenient to write the partonic cross- 
section in terms of a dimensionless coefficient function 
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FIG. 1: The 0(a 3 ) contribution to the coefficient function 
Eq. Q as a function of N, for tuh = 125 GeV (upper plot) 
and ran = 600 GeV (lower plot). In each case we show the ex- 
act result and the pointlike and logarithmic approximations, 
as well as the ratio of the latter two to the exact result. 



C(z,a s ) 



a(z,a s ) = cr zC(z,a s ), 



(3) 



where cto is the LO partonic cross-section, so the gg co- 
efficient function has an expansion in powers of a s 

C(z,a s ) = 5(l-z) + a s C {1) (z) + a 2 s CW(z) + 0(a 3 s ). (4) 

Upon Mellin transformation 

JV-l ffawj?) 



a(N,m 2 H ) = / dTT f 

Jo T 

C(N,a s )= [ dzz N - 1 C(z,a s ) 7 



dzz N ~ x ££(z,y?) 



(5) 
(6) 
(7) 



the convolution in Eq. ([I]) turns into an ordinary product: 

a(N, m% ) = (7 J?(N, m 2 H ) C(N,a s ). (8) 

The Mellin transformation maps the large r — > 1 region 
into the large N — > oo region, and the small r — > region 
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FIG. 2: Same as Fig.[TJ but to 0(a 2 s 



into the small N —> N s region, with N s the rightmost sin- 
gularity of ^(A^, m 2 H ) (i.e. the convergence abscissa of the 
Mellin transform). For gluon-initiated processes, N s = 1 
to all perturbative orders. 

The dominant partonic kinematic region can then be 
determined through a saddle point argument, by com- 
puting the value of N which provides the dominant con- 
tribution to the Mellin inversion integral: 



ct(t, m%) 



1 



c+zoo 



dNr- N a(N, m 2 H ) 



Namely, we define 



(c > N s ). 
(9) 



E(N,T,m 2 H ) =N In- + In a(N, m 2 H ), (10) 
r 



so that 



ct(t, m 2 H ) 



1 

27ri 



dNe E{N,r, m l)_ (n) 



In the saddle point approximation, the integral Eq. 

is dominated by the value of the exponent E(N,T,mfP) 

at its stationary point iVo 



dE(N,T,m 2 H ) 



dN 



= 0, 



(12) 



N=N n 
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and by the behaviour in its vicinity. The position of the 
saddle point is a function No = /Vo(t, m^), solution of 
Eq. ( 12 ). Hence, for any value of the physical kinematics, 



the questions whether the partonic cross-section is well 
approximated by its pointlike limit or resummation is 
advantageous are answered by verifying whether this is 
the case for C(N , a s ). 

A unique real saddle point is present due to the drop 
of the cross-section a(N,m 2 H ) as N grows. This drop 
is driven by the parton luminosity J£(N, /i 2 ), which thus 
controls the position of the saddle N . The drop of the lu- 
minosity at large TV (and its growth at small N) reflects in 
turn the drop of the PDFs and luminosity as z — > 1 (and 
their growth as z — > 0). However, for large N, C(N,a s ) 
actually grows with N. This growth, which is due to 
the logarithmically enhanced contributions, is only pos- 
sible because the partonic cross-section is a distribution, 
rather than an ordinary function: the Mellin transform 
of an ordinary positive function is easily proven to be a 
decreasing function of N. Nevertheless, the parton lumi- 
nosity always offsets this increase, because the physical 
cross-section <t(t, m 2 H ) is an ordinary positive function 
and thus must decrease with N. The faster the growth of 
the cross-section, the better the saddle point approxima- 
tion, which is thus especially good for gluon-dominated 
processes, as the gluon is more peaked at small z than the 
quark, and thus the gluon-gluon luminosity drops faster 
than quark luminosity. 

The stationary point is determined by the interplay 
in Eq. ( 10 ) of the rise of the first term and the drop 



of the hadronic cross-section a(N,rnjj), which in turn is 
the product of the coefficient function and the luminosity. 
The shape of the NLO and NNLO contributions to the co- 
efficient function is shown in Figs. [T] and [2j for two values 
of the Higgs mass which are allowed by present data [2J. 
Clearly, when the partonic cross-section rises monoton- 
ically with N, it is the drop of the parton luminosity 
which determines the drop of the hadronic cross-section 
and thus the position of the saddle, and even when it 
drops, the decrease of its hadronic counterpart is much 
stronger in the presence of a luminosity, so the location 
of the saddle is substantially larger. 

We have determined the position of the saddle at NLO 
and NNLO, using NNPDF2.1 parton distributions [5] at 
the corresponding order. Results are shown in Fig. [3j 
which is the main result of this paper. The position 
of the saddle N depends on two independent variables, 
which can be chosen out of the three variables mn-, s, 
and r = m 2 H / s. If results for N are shown as a function 
of r, the dependence on the other variable (be it mji or 
s) becomes very slight, because it enters only through the 
scale dependence of a s and the parton distributions: this 
is explicitly seen in the top plot of Fig. [3j where results 
are shown for two different values of m#; the dependence 
would be similar if y/s were varied instead by a compara- 
ble factor. The dependence of the position of the saddle 
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FIG. 3: The position of the saddle point No for the Mellin 
inversion integral Eq. (J9j) as a function of: r = m 2 H /s (top); 
win for three different values of the collider energy s (middle) ; 
s for two different values of mn (bottom). The curves in 
the top plot depend very weakly on either ttih (shown) or 
s. Both the NLO and NNLO are shown, computed using the 
appropriate NNPDF2.1 PDF set [5] in each case. 



on the perturbative order is completely negligible (as also 
shown), and so is the dependence on whether the point- 
like approximation is used or not. We have also checked 
that varying the rcnormalization and factorization scales 
by a factor two changes the value of No by 1% or less, 
even though more dramatic scheme changes which affect 
the infrared behaviour of the PDF such as suggested in 
Ref. [Hj might have somewhat larger effects. The impact 
of changing the PDF set or the value of a s in a reasonable 
range is rather less than that. 

The size of the region around the value of No which 
dominates the integral can be estimated by computing 
the second derivative of \na(N,m'j I ), which gives the 
width of the gaussian integral which approximates the 
Mellin inversion in the complex N plane. We find that 
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a one-sigma region corresponds to a variation of Nq by 
about 25% about its central value. So, within the accu- 
racy of our saddle point approximation, the distinction 
between different curves at fixed r is of little import. For 
clarity and completeness, in Fig. [3] we show the position 
of the saddle as a function of the Higgs mass for the three 
values of s relevant for the LHC, and as a function of the 
center of mass energy for two values of the Higgs mass. 

We can now assess both the adequacy of the pointlike 
approximation, and the desirability of resummation for 
given values of s and to#: first, using Fig. [3] the given 
hadronic kinematics can be translated into a value of 
Ao. Then, for this value of Nq we can check whether 
the pointlike approximation is accurate, and threshold 
resummation is advisable. To this purpose, in Figs. [T]- 
[2] we compare the exact NLO and NNLO contributions 
to the coefficient function to its various approximations. 
At NLO, we use the implementation of Ref. [TU] of the 
original exact result of Refs. [TT1 H2; at NNLO a full 
exact computation is not available, so for light Higgs we 
use the expansion of Ref. [4] matched to the exact small 
z limit of Ref. [T3], while for heavy Higgs, where the 
expansion of Ref. [I] is unstable, the "finite mt" curve of 
Fig. [2] is merely given by correcting the pointlike result 
through the inclusion of the first order correction in -^ it - 
from Ref. |4j (and it thus only indicates the location of the 
region where the finite m t corrections are likely relevant). 

The pointlike results at NLO O [15] and NNLO pj] 
have been long known. As a soft part, we show all con- 
tributions which survive the A~ — > 00 limit, defined as 
the exact Mellin transform of all contributions of the 

form z — J and all contributions proportional 

to 5(1 — z) to the NLO and NNLO coefficient functions. 
While the coefficients of these terms are fixed by soft re- 
summation, there is a certain latitude in defining which 
subleading terms to include. Our choice reproduces the 
exact soft kinematics [7j[17] thereby optimizing the agree- 
ment with the exact result to all orders. The region in 
which the soft limit defined in this way is close to the full 
result is the region in which one expects soft resumma- 
tion to improve the accuracy of the computation, even 
when all-order resummation is not mandatory. 

Inspection of Figs. [T[j2] shows that the region in which 
the pointlike approximation is adequate is very close to 
the region in which the soft approximation is good, so 
indeed the success of the former might be explained by 
the accuracy of the latter. Both at NLO and NNLO the 
relevant region is roughly N > 2. This is the region 
where the partonic cross-section starts to rise with N, 
driven by logarithmically enhanced contributions. The 
apparent failure of both approximations for heavy Higgs 
at NNLO is likely to be due to the fact that in this case 
the "finite m t " computation is in fact only approximate. 

The fact that the same behaviour is observed at NLO 
and NNLO is not accidental. On the one hand, the po- 



sition of the saddle is largely determined by the parton 
luminosity, which is perturbatively very stable at NLO 
and beyond [8|. On the other hand, the shape of the 
coefficient functions and the dominance of soft terms are 
mostly controlled by the location of the leading small- 
and large- A~ singularities, which can be checked to be sta- 
ble even upon all-order resummation 7J. This supports 
the expectation that the desirability of resummation and 
the all-order reliability of the pointlike approximation can 
be assessed on the basis of the known low orders. 

At the LHC, r is quite small: if ran ~ 125 GeV, 
r - 1CT 4 and if m H ~ 600 GeV, r - 10~ 3 - 10~ 2 , 
thereby leading to values of Ao, close to the transition 
value N ^ 2, for which resummation is at best desirable, 
but certainly not mandatory, as a s In Ao < 1. Never- 
theless, Fig. [3j in which the N = 2 line has been drawn 
for ease of reference, shows that for a heavy Higgs with 
rriH 600 GeV the pointlike approximation is adequate, 
and Sudakov resummation clearly advantageous, for any 
collider with center of mass energy up to y/s < 30 TeV, 
and thus certainly at the LHC. 

Even for a light Higgs with m# ~ 125 GeV the point- 
like approximation is fine for a collider with energy of 
7 TeV, but it may start failing as the energy is raised, 
and becomes inadequate at the LHC with y/s — 14 TeV, 
where the saddle drops at Vo ~ 1.9 so the finite-mass 
corrections to Fig. [TJreach the percent level, though 
care should be taken because, close to the region where 
the approximation breaks down, the conclusion may de- 
pend on various details and approximations. However, 
because the breakdown of the pointlike approximation is 
in significant part due to the presence of spurious high- 
energy double logs [15] . it is expected to be more notice- 
able in less inclusive observables. Correspondingly, with 
this higher center of mass energy, Sudakov resummation 
is likely to stop being advantageous for a light Higgs. 

In summary, we have provided a way of assessing 
whether the pointlike approximation is adequate and 
whether Sudakov resummation is desirable based on the 
behaviour of the known first few orders of the partonic 
cross-section. Based on it, we conclude that for a light 
Higgs with mji ~ 125 GeV the full inclusion of finite 
top mass corrections to its cross-section is likely to be 
important for accurate phenomenology at the LHC with 
= 14 TeV. 

We thank F. Tackmann for a critical reading of the 
manuscript, and R. Harlander and A. Vicini for providing 
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